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We have derived consistent sets of band parameters (band gaps, crystal field-splittings, band 
gap deformation potentials, effective masses, Luttinger and i?p parameters) for AIN, GaN, and 
InN in the zinc-blende and wurtzite phases employing many-body perturbation theory in the 
GqWo approximation. The GqWo method has been combined with density-functional theory (DFT) 
calculations in the exact-exchange optimized effective potential approach (OEPx) to overcome the 
limitations of local-density or gradient-corrected DFT functionals (LDA and GGA). The band struc- 
tures in the vicinity of the F-point have been used to directly parameterize a 4x4 k-p Hamiltonian 
to capture non-parabolicities in the conduction bands and the more complex valence-band structure 
of the wurtzite phases. We demonstrate that the band parameters derived in this fashion are in 
very good agreement with the available experimental data and provide reliable predictions for all 
parameters which have not been determined experimentally so far. 
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I. INTRODUCTION 

The group Ill-nitridcs AIN, GaN, and InN and their 
alloys have become an important and versatile class of 
semiconductor materials, in particular for use in opto- 
electronic devices and high-power microwave transistors. 
Current applications in solid state lighting [light emitting 
diodes (LEDs) and laser diodes (LDs)] range from the 
visible spectrumii^i^ to the deep ultra-violet (UV}^i^, 
while future applications as, e.g., chemical sensorsii^i^ or 
in quantum cryptographj*iS are being explored. 

For future progress in these research fields reliable 
material parameters beyond the fundamental band gap, 
like effective electron masses and valence-band (Lut- 
tinger or Luttinger-like) parameters, are needed to aid 
interpretation of experimental observations and to simu- 
late (hetero-)structures, like, e.g., optoelectronic devices. 
Material parameters can be derived from first-principles 
electronic-structure methods for bulk phases, but the size 
and complexity of structures required for device simula- 
tions currently exceeds the capabilities of first-principles 
electronic-structure tools by far. To bridge this gap first- 
principle calculations can be used to parameterize sim- 
plified methods, like the k-p method j^^i^^'^'^i^'^ the empir- 
ical tight-binding (ETB) method ^^li^iiiii-^ or the empir- 
ical pseudo-potential method (EPM)i^, which are appli- 
cable to large-scale heterostructures at reasonable com- 
putational expense. 

In this Article we use many-body perturbation theory 
in the GqWo approximation,^*^ — currently the method of 
choice for the description of quasiparticlc band struc- 
tures in solids^iiS^i^^ — in combination with the k • p 
approachfiiii^ii^iii^ to derive a consistent set of material 
parameters for the group-Ill nitride system. The k • p - 



Hamiltonian is parameterized to reproduce the GqWq 
band structure in the vicinity of the F-point. Since the 
parameters of the k • p method are closely related or, 
in some cases, even identical to basic band parameters, 
many key band parameters can be directly obtained us- 
ing this approach. 

The k-p -model Hamiltonian is typically parameterized 
for bulk structures and is then applicable to heterostruc- 
tures with finite size (e.g., micro- and nanostructures) 
within the envelope-function scheme.^"* Ideally, the pa- 
rameters are determined entirely from consistent experi- 
mental input. For the group-III-nitrides, however, many 
of the key band parameters have not been conclusively 
determined until now, despite the extensive research ef- 
fort in this fieldi^iSS. In a comprehensive review Vur- 
gaftman and Meyer summarized the field of III-V semi- 
conductors in 2001 and recommended up-to-date band 
parameters for all common compounds and their alloys 
including the nitridesi^i Only two years later they real- 
ized that it is striking how many of the nitride proper- 
ties have already been superseded, not only quantitatively 
but qualitatively!^ They proceeded to remedy that obso- 
lescence, by providing a completely revised and updated 
description of the band parameters for nitride- containing 
semiconductors in 2003.^^ While this update includes ev- 
idence supporting a revision of the band gap of InN from 
its former value of 1.9 eV to a significantly lower value 
around .7 eV r^^^^-i^^^^^ they had to concede that in 
many cases experimental information on certain param- 
eters was simply not availablci^ This was mostly due 
to growth-related difficulties in producing high quality 
samples for unambiguous characterization. In the mean- 
time the quality of, e.g., wurtzite InN samples has greatly 
improved^^ and even the growth of the zinc-blende phase 
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has advanced Nevertheless, many of the basic material 
properties of the group-III nitrides are still undetermined 
or, at least, controversial. 

On the theoretical side, certain limitations of density- 
functional theory (DFT) in the local-density or gen- 
eralized gradient approximation (LDA and GGA, 
respectively) — currently the most wide-spread ab-initio 
electronic-structure method for poly-atomic systems — 
have hindered an unambiguous completion of the miss- 
ing data. To overcome these deficiencies we use 
Go Wq calculations based on DFT calculations in the 
exact-exchange optimized effective potential approach 
(OEPx) to determine the basic band parameters. We 
have previously shown that the OEPx-|-GoWo approach 
provides an accurate description of the quasiparticle band 
structure for GaN, InN and II- VI compounds^i^i^^i^. 
The quasiparticle band structure in the vicinity of the 
F-point is then used to parameterize a 4 x 4 k-p Haniilto- 
nian to determine band-dispersion parameters, like effec- 
tive masses, Luttinger parameters, etc. This allows us to 
take the non-parabolicity of the conduction band, which 
is particularly pronounced in In N^^'^"^ , and the more com- 
plex valence band structure of the wurtzite phases into 
account properly. 

This paper is organized as follows: In Section |ll] we 
briefly introduce the GqWq approach and its application 
to the group-III nitrides, followed by a discussion of cer- 
tain key parameters of the quasiparticle band structure, 
such as the fundamental band gaps (and their depen- 
dence on the unit-cell volume) and the crystal-field split- 
ting energies (Section IIIip . In Section |IV] we present 
our recommendations for the band dispersion parame- 
ters (k • p parameters) of the wurtzite and zinc-blende 
phases of AIN, GaN and InN. A detailed discussion of 
the parameter sets is given in Section [IVB] together with 
a comparison to experimental values and parameter sets 
obtained by other theoretical approaches. Our conclu- 
sions are given in Sec. |Vl 



II. QUASIPARTICLE ENERGY 
CALCULATIONS 



A. GW based on exact-exchange DFT 

The root of the deficiencies in LDA and GGA for de- 
scribing spectroscopic properties like the quasiparticle 
band structure can be found in a combination of different 
factors. LDA and GGA are approximate (jellium-based) 
exchange-correlation functionals, which suffer from in- 
complete cancellation of artificial self- interaction and lack 
the discontinuity of the exchange-correlation potential 
with respect to the number of electrons. As a con- 
sequence the Kohn-Sham (KS) single-particle eigenval- 
ues cannot be rigorously interpreted as the quasiparticle 
band structure as measured by direct and inverse pho- 
toemission. This becomes most apparent for the band 
gap, which is severely underestimated by the Kohn-Sham 




FIG. 1: LDA Kohn-Sham calculations incorrectly predict 
wurtzite InN to be a metal with the wrong band ordering 
at the r point. In OEPx(cLDA) the band gap opens and InN 
correctly becomes a semiconductor, thus providing a more 
suitable starting point for subsequent quasiparticle energy cal- 
culations in the GqWo approximation. 



eigenvalue difference in LDA and GGA. For InN this even 
results in an overlap between the conduction and the va- 
lence bands and thus an effectively metallic state, as dis- 
played in Fig. [H It goes without mentioning that a k-p 
parameterization derived from this LDA band structure 
would not appropriately reflect the properties of bulk 
InN. 

Many-body perturbation theory in the GW approaches 
presents a quasiparticle theory that overcomes the defi- 
ciencies of LDA and GGA and provides a suitable de- 
scription of the quasiparticle band structure of weakly 
correlated solids, like AIN, GaN and lnN^L22^ Most 
commonly, the Green's function Go and the screened po- 
tential Wo required in the GW approach (henceforth de- 
noted GqWo) are calculated from a set of DFT Kohn- 
Sham single particle energies and wave functions. The 
DFT ground state calculation is typically carried out in 
the LDA or GGA and the quasiparticle corrections to 
the Kohn-Sham eigenvalues are calculated in first order 
perturbation theory (LDA/GGA -t- GqWq) without re- 
sorting to self-consistency in G and Wi^ 

While the LDA-l-GoWo approach is now almost rou- 
tinely applied to bulk materials r^i^SiS^ GqWq calcu- 
lations for GaN and InN have been hampered by 
the deficiencies of the LDA. For zinc-blende GaN the 
LDA+GqWo band gap of 2.88 eV^^ is still too low 
compared to the experimental 3.3 eV^^i^i^ while for 
InN the LDA predicts a metallic ground state with incor- 
rect band ordering. A single Go Wo iteration proves not 
to be sufficient to restore a proper semiconducting state 
and only opens the band gap to 0.02 - 0.05 cV,^^'"*^ which 
is still far from the experimental value of ~0.7eV.— i^Sii^ 

Here we apply the GqWo approach to DFT calcu- 
lations in the exact-exchange optimized effective po- 
tential approach (OEPx or OEPx(cLDA) if LDA cor- 
relation is included). In contrast to LDA and GGA 
the OEPx approach is fully self-interaction free and 
correctly predicts InN to be semiconducting with the 
right band ordering in the wurtzite phaseS^^ as Fig. [1] 
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ao [A] 


CO [A] 


C()/ao 


u 


zb-AlN 


4.370 








zb-GaN 


4.500 








zb-InN 


4.980 








wz-AlN 


3.110 


4.980 


1.6013 


0.382 


wz-GaN 


3.190 


5.189 


1.6266 


0.377 


wz-InN 


3.540 


5.706 


1.6120 


0.380 



TABLE L Experimental lattice parameters (ao, co and u) 
adopted in this work (see text), for zinc blende (zb) and 
wurtzite (wz) AIN, GaN, and InN.— 



demonstrates. For II- VI compounds, GaN and ScN we 
have previously illustrated the mechanism behind this 
band gap opening in the OEPx approach)22iMi2^i^ii2 
which brings the OEPx Kohn-Sham band gap much 
closer to the experimental one. Combining OEPx with 
GqWq^ in this fashion yields band gaps for II- VI com- 
pounds, Ge, GaN and ScN in very good agreement with 
experimenti^i^ii^^i^^dl For wurtzite InN the band gap of 
0.7 eV and the non-parabolicity of the conduction band 
(CB)^ (shown in Fig. [1]) strongly supports the recent 
experimental findings.— '^'^i^'^ In addition we have 
shown that the source for the startling, wide interval of 
experimentally observed band gaps can be consistently 
explained by the Burstein-Moss effect (apparent band 
gap increase with increasing electron concentration in the 
conduction band) by extending our calculations to finite 
carrier concentrations.— 



B. Computational Parameters 

The LDA and OEPx calculations in the present work 
were performed with the plane-wave, pseudopotential 
code S/PHI/nX,— while for the GqWq calculations we 
have employed the GqWo space-time method^^ in the 
gwst implementationi^i^i^^ Local LDA correlation is 
added in all OEPx calculations. Here we follow the 
parametrization of Perdew and Zunger— for the corre- 
lation energy density of the homogeneous electron gas 
based on the data of Ceperley and Alder This com- 
bination will in the following be denoted OEPx(cLDA). 
Consistent pseudopotentials were used throughout, i.e. 
exact-exchange pseudopotentials^^ for the OEPx(cLDA) 
and LDA ones for the LDA calculations. The cation 
d-electrons were included explicitly.^-"'*- For additional 
technical details and convergence parameters we refer to 
previous work.— 



C. Lattice Parameters 

All calculations are carried out at the experimental 
lattice constants reported in Tab. [J and not at the ab 
initio ones to avoid artificial strain effects in the derived 
parameter sets. For an ab initio determination of the 



lattic constants consistent with the GqWo calculations 
the crystal structure would have to be optimized within 
the Go Wo formalism, too. However, Go W^o total en- 
ergy calculations for realistic systems have, to our kon- 
wledge, not been performed, yet, and the quality of the 
Go Wo total energy for bulk semiconductors has not been 
assessed so far.^?- The alternative ab initio choices, LDA 
and OEPx(cLDA), give different lattice constants^ and 
would thus introduce uncontrolable variations in the cal- 
culated band parameters that would aggrevate a direct 
comparison. 

The thermodynamically stable phase of InN at the 
usual growth conditions is the wurtzite phase. Reports 
of a successful growth of the zinc-blende phase have been 
scarce. Recently, high-quality films of zb-InN grown on 
indium oxide have been obtained by Lozano et al^ We 
adopt their lattice constant of 4.98 which is in good 
agreement with previous reports of 4.98 4.986 A^— 
and 5.04 AS for wz-InN grown on different substrates. 
For zb-AIN reports of successful growth are even scarcer. 
Petrov et al. first achieved to grow AIN in the zinc-blende 
phase and reported a lattice constant of ao=4.38 Ai^ 
This was later refined by Thompson et al. to a=4.37 
which is the value we adopt in this work. For zb-GaN we 
follow the work of Lei et al^^ and chose ao=4.50 A. 

Although wurtzite is the phase predominantly grown 
for InN, reported values for the structural parameters 
still scatter appreciablyi^ In order to determine the 
effect of the lattice constants on the band gap {Eg) 
and the crystal-field splitting (Acr) we have explored 
the range between the maximum and minimum val- 
ues of ao and cq/qq reported in Ref. [6^ by perform- 
ing OEPx(cLDA)+GoWo calculations at the values listed 
in Tab. [iTl Since u remains undetermined in Ref. [6^ 
we have optimized it in the LDA. Neither u, i?g, nor 
Acr depends sensitively on the lattice constants in this 
regime and we have therefore adopted the mean values of 
ao=3.54 A, co=5.706 A(co/ao=1.612) and u=0.380 (the 
LDA-optimized value) for the remainder of this article. 

For WZ-A1N and wz-GaN the lattice constants are more 
establishedi^iS For wz-AIN we adopt Schulz and Thier- 
mann's values of a=3.110 A, c=4.980 A, and u=QM2^ 
which are close to those reported by Yim et al^ Schulz 
and Thiermann also provide a value for the internal pa- 
rameter u, which is identical to the one we obtain by 
relaxing u in the LDA at the experimental ao and co 
parameters. The same is true for wz-GaN. Schulz and 
Thiermann's values of a=3.190 Aand c=5.189 A™ are 
close to those first reported by Maruska and Tietjenj^^ 
but in addition offer a value of M=0.377, which corre- 
sponds to our LDA-relaxed value at the same lattice pa- 
rameters. Note, that the lattice parameters of wz-InN 
and wz-GaN have been refined compared to our recently 
published calculations.'^'* The influence of the adjustment 
on the different band parameters will be discussed where 
necessary. 
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ao (A) 






OEPx(cLDA)+G'oWo 


OEPx(cLDA) 


co/ao 


u 


AcR (eV) 


E, (eV) 


AcH, (eV) 


E, (eV) 


3.535 


1.612 


0.380 


0.067 


0.71 


0.079 


1.01 


3.540 


1.612 


0.380 


0.066 


0.69 


0.079 


1.00 


3.545 


1.612 


0.380 


0.065 


0.68 


0.079 


0.98 


3.540 


1.611 


0.380 


0.065 


0.70 


0.078 


1.00 


3.540 


1.612 


0.380 


0.066 


0.69 


0.079 


1.00 


3.540 


1.613 


0.380 


0.068 


0.69 


0.081 


0.99 



TABLE II: Band gap (-Eg), and crystal field splitting (Acr) for wurtzite InN in tiie range of tiie experimentally reported values 
of the structural parameters ao and co/ao (Ref. i68i' ) and u determined in LDA for ao=3.54 Aand co/ao=1.612. 



III. BAND GAPS, CRYSTAL-FIELD 
SPLITTINGS, AND BAND-GAP DEFORMATION 
POTENTIALS 

We will now discuss the quasiparticle band structure 
of AIN, GaN and InN in their zinc-blende and wurtzite 
phases in terms of certain key band parameters such 
as the band gap (-Eg), the crystal field spUtting (Acr) 
in the wurtzite phase and the band-gap volume defor- 
mation potentials ay- At the end of this section we 
will draw a comparison between LDA and OEPx(cLDA) 
based GqWq calculations for AIN. 



A. Band Gaps 

The OEPx(cLDA)+GoPVo band gaps for the three ma- 
terials and two phases are reported in Tab. IIIII together 
with the LDA and OEPx(cLDA) values for compari- 
son. For GaN and InN the OEPx(cLDA)-hGo Wo band 
gaps have been reported previously in Ref. [3J. There 
we have also argued that the wide interval of exper- 
imentally observed band gaps for InN can be con- 
sistently explained by the Burstein-Moss effect. The 
OEPx(cLDA)-HGo Wo value of 0.69 eV for wz-InN^s sup- 
ports recent observations of a band gap at the lower end 
of the experimentally reported range. For zinc-blende 
InN, which has been explored far less experimentally, 
our calculated band gap of 0.53 eV also agrees very 
well with the recently measured (and Burstein-Moss cor- 
rected) 0.6 eV;^ 

For GaN the band gaps of both phases 
are well established experimentally and our 
OEPx(cLDA)-|-Go Wo calculated values of 3.24 eVM 
and 3.07 eV agree to within 0.3 eV. 

For AIN experimental results for the band gap of 
the wurtzite phase scatter appreciable, whereas for zinc 
blende only one value has - to the best of our knowl- 
edge - been reported so far. Contrary to GaN, the 
OEPx(cLDA)-|-GoWogaps for AIN are larger than the 
experimentally reported values. 



B. Crystal-Field Splitting 



Experimental values for the crystal-field splitting, 
Acr, of wz-GaN scatter between 0.009 and 0.038 eV 
(Tab. Unj. The OEPx(cLDA)-hGo Wo value of 0.033eV 
supports a crystal-field splitting within this range. 

Theoretical^ and experimental^^ investigations of wz- 
AIN agree upon the fact that the crystal-field splitting of 
AIN is negative. Our calculations also yield a negative 
value of Acr — —0.295 eV. This result supports a crystal- 
field splitting in AIN below — 0.2eV, as reported by 
Chen et ai.,S, rather than a small negative value between 
—0.01 and —0.02, as implied by the results of Freitas et 
al^ For wz-InN a crystal-field splitting between 0.019 
and 0.024 eV has been reported recently^^ This value is 
significantly smaller than the OEPx(cLDA)-|-GoWo value 
of 0.07 eV. 

The crystal-field splitting is known to be sensitive to 
lattice deformations, such as changes in the co/ag ratio 
or the internal lattice parameter ^, ,89,90,91 Therefore, the 
discrepancy between experiment and theory might stem 
from the uncertainties of the lattice parameters of wz- 
InN (cf Sec. Ill C|) . However, varying the co/ao ratio or 
the unit cell volume within the experimental range dis- 
cussed in Section III CI yields values for Acr which are 
always larger than 0.06 eV (Tab. |TT|), leaving only the 
internal lattice parameter u as possible source of error. 
This parameter is - at least for GaN - known to have a 
large influence on the crystal-field splittingi^ Although 
the LDA-optimized u values are in very good agreement 
with experimental values for GaN and AIN, experimental 
confirmation of the u parameter of InN is still pending. 
We therefore calculated the crystal-field splitting of wz- 
InN for different values of u (and ao and cq fixed at the 
values listed in Tab. [l| between 0.377 and 0.383. Gen- 
erally, Acr decreases with increasing u, but even for u 
as large as 0.383, the crystal-field splitting is still larger 
than 0.05 eV The discrepancy between the experimen- 
tal report and the OEPx(cLDA)-|-GoWo calculations can 
hence not be attributed to the uncertainties in the lat- 
tice parameters and has to remain unsettled for the time 
being. 
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par am. 


OEPx(cLDA)+GoW"o 
(this work) 


Exp. 


LDA 
(this 


OEPx(cLDA) 
work, for comparison) 


wurtzitc 


wz-AlN 




6.47 


6.0-6.3" 


4.29 


5.73 




AcR 


-0.295 


-0.230" 


-0.225 


-0.334 


wz-GaN 




3.24 


3.5" 


1.78 


3.15 




AcR 


0.034 


0.009-0.038" 


0.049 


0.002 


wz-InN 




0.69 


0.65-0.8^ 




1.00 




AcR 


0.066 


0.019-0.024'' 




0.079 


zinc blende 


zb-AlN 




6.53 




4.29 


5.77 






5.63 


5.34" 


3.28 


5.09 


zb-GaN 




3.07 


3.3^* 


1.64 


2.88 


zb-InN 


Sg 


0.53 


0.6^ 




0.81 



"Reference 1 65l 

''References VfW^ [tI [TsI. [T^. [TTI. andlzi 

'^Reference |73| . 
■^Reference s [ill, [il and [H 

''Reference |26| and references therein. 
■^Reference 1791 

^References lii, H Ell, S IH, [H, and IH 

Reference [sl 

TABLE III: Band gaps {Eg) and crystal-field splittings (Acr)^ for the wurtzite and zinc-blende phases of AIN, GaN, and InN. 
All values are given in eV. 



OEPx(cLDA)+GoW^o LDA OEPx(cLDA) 







(this work) 


(this 


work, for compar. 


wurtzitc 


AIN 


wz 
Oly 


-9.8 


-8.8 


-8.9 


GaN 


„ wz 

Q-V 


-7.6 


-6.8 


-6.5 


InN 




-4.2 




-3.0 


zinc blende 


AIN 


zb I'-l' 


-10.0 


-9.1 


-9.1 






-1.8 


-0.6 


-0.9 


GaN 


at 


-7.3 


-6.4 


-6.1 


InN 




-3.8 




-2.6 



TABLE IV: Volume band gap deformation potentials (av) for 
the wurtzite and zinc-blende phases of AIN, GaN, and InN. 
The volume deformation potentials can be transformed into 
pressure deformation potentials using the bulk moduli of the 
respective materials (see text). All values are given in eV. 



C. Band Gap Deformation Potentials 

For the hydrostatic band gap deformation potentials 
the band gaps have been calculated at different vol- 
umes {V) between ±2 % around the equilibrium vol- 
ume Vq. In the explored volume range the band gaps 
vary linearly with h\{V/Vo)- The linear coefficient is 
then taken as the hydrostatic volume deformation po- 
tential ay- The calculated band gap deformation po- 
tentials are listed in Tab. |TV] for LDA, OEPx(cLDA), 
and OEPx(cLDA)+GoM^o ■ We observe that for all com- 
pounds and phases the quasiparticle deformation po- 
tential is larger in magnitude than that of the DFT 



(LDA and OEPx(cLDA)) calculations, i.e., the band 
gaps vary stronger with volume deformations. Hydro- 
static band gap deformation potentials obtained from 
LDA+U calculations have recently been reported for the 
wurtzite phases of GaN and InN^^. Unlike in OEPx, 
where an improved description of the p-d hybridization 
is achieved by the full removal of the self-interaction 
for all valence states, the LDA-I-U approach reduces 
the p-d repulsion by adding an on-site Coulomb corre- 
lation U only to the semicore c?-electrons. With refer- 
ence to the OEPx(cLDA)-f GoWo deformation potentials, 
the LDA+U improves upon LDA for GaN (a^°^+u=- 
7.7 eV), but worsens for InN ^ 
4.2 eV in ReL W 



=-3.1 eV, ai>^^= 



Experimentally the band gap deformation potential 
is usually measured as a function of the applied pres- 
sure, which aggravates a direct comparison to our cal- 
culated volume deformation potentials. However, since 
B = —dP/d\nV, where B is the bulk modulus and P 
the pressure, the pressure deformation potential ap can 
be expressed in terms of ay according to ap = —ay/B. 

Experimentally reported values for the bulk modulus of 
wz-GaN scatter between 1880 and 2450 kbar,23i94^j96^ 
Using these values, our volume deformation potential 
of av'=-7.6 eV would translates into a pressure defor- 
mation potential in the range of 3.1 - 4.0 meV/kbar, 
which is comparable to the experimentally determined 
range of 3.7 and 4.7 meV/kbar,2L98i99iioOiioiii02 r^j^-^ 

large uncertainty has been partially ascribed to the low 
quality of earlier samples and substrate-induced strain 
effects.i^ The fact that the pressure dependence of the 
band gap is sublinear (unlike the volume dependence) 
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further questions the accuracy of hnear or quadratic fits 
for the extraction of the deformation potentials in the 
experiments 

For wz-InN experimentally reported values are sparse. 
Franssen et al. determined a hydrostatic pressure de- 
formation potential of 2.2 meV/kbar-i2^, while Li et al. 
found 3.0 mcV/kbarii^ This range agrees with our the- 
oretical one of 2.8 - 3.3 meV/kbar, using for the con- 
version of volume to pressure deformation potentials the 
bulk modulus range of 1260 - 1480 kbar— quoted in 
the literature. 

For wz-AlN we are only aware of one experi- 
mental study reporting a pressure deformation po- 
tential of 4.9 meV/kbari^ With experimental bulk 
moduh between 1850 and 2079 kbari^sooeaoi ^j^g 
OEPx(cLDA)-|-Go Wo pressure deformation potential of 
WZ-A1N would fall between 4.7 and 5.3 meV/kbar strad- 
dling the experimentally reported value. 

To our knowledge, no experimental information on the 
deformation potential of zb-AlN and zb-InN are avail- 
able. For zb-GaN our computed volume deformation 
potential of ay=-7.3 eV translates to a pressure defor- 
mation potential range of 3.0 - 3.9 meV/kbar using the 
same bulk modulus range as for wz-GaN. This range is 
slightly below the experimentally reported range of 4.0 
- 4.6 meV/kbarii22ii2^ii2S Employing a semi-empirical 
approach to overcome the band gap underestimation of 
the LDA (LDA-plus-correction (LDA-I-C), see also dis- 
cussion in Section IIVB ip , Wei and Zunger found vol- 
ume deformation potentials of — 10.2eV [zb-AIN (F — 
F)], -l.leV [zb-AIN {T - X)], -7.4eV (zb-GaN) and 
— 3.7eV (zb-InN)^^" in good agreement with our full 
OEPx(cLDA)-HGo Wo calculations (see Tab.HVl). 



D. Comparison between LDA+Go Wo and 
OEPx(cLDA)+ Go Wo 

For the materials presented in this article a meaning- 
ful comparison between LDA and OEPx(cLDA) based 
Go Wo calculations can only be constructed for AIN for 
reasons given in Section Hi Al Figure [2] displays the band 
structure of wz-AIN in the four approaches discussed in 
this Article. The "band gap problem" has been elimi- 
nated from this comparison by aligning the conduction 
bands at the minimum of the lowest conduction band 
(ecbm) and the valence bands at the maximum of the 
hightest valence band (evBivi)- For this large gap ma- 
terial three main conclusions can be drawn from Fig. [21 
First, LDA and both Go Wo calculations yield very sim- 
ilar band dispersions. Or in other words the GqWo cor- 
rections to the LDA in the LDA+GqWo approach are 
not k-point dependent shifting bands almost rigidly. A 
rigid shift between conduction and valence bands is fre- 
quently referred to as "scissor operator" . Fig.[2l however, 
illustrates that this shift is not identical for all bands, 
which cannot be attributed to a single scissor operator. 
Second, the dispersion obtained in OEPx(cLDA) devi- 
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FIG. 2: Comparison between LDA (solid lines), 
LDA -f Go Wo (open circles), OEPx(cLDA) (dashed lines) 
and OEPx(cLDA)-(-GoWo (solid circles) for wz-AlN. In both 
GoWo calculations bands are essentially shifted rigidly 
compared to the LDA, whereas OEPx(cLDA) yields bands 
with different dispersion. 



ates from the other three approaches, which is consistent 
with the observation made for wz-InN in Fig. [TJ We at- 
tribute this behaviour to the approximate treatment of 
correlation in the OEPx(cLDA) approach and the fact 
that the band structure in OEPx(cLDA) is a Kohn- 
Sham and not a quasiparticle band structure. While 
the LDA benefits from a fortuitous error cancellation be- 
tween the exchange and the correlation par1>iii, this is no 
longer the case once exchange is treated exactly in the 
OEPx(cLDA) scheme. Using a quasiparticle approach 
with a more sophisticated description of correlation, like 
the GW method, then notably changes the dispersion 
of the OEPx(cLDA) bands. As we will demonstrate 
in the next Section this will lead to markedly different 
band parameters not only for the conduction but also 
for the valence bands (cf Tab. IVII|) . Against common 
believe OEPx(cLDA) calculations without subsequent 
Go Wo calculations may therefore provide a distorted pic- 
ture and we would advise against deriving band parame- 
ters from OEPx or OEPx(cLDA) band structures alone. 
Third, unlike in the LDA-I-GqWo case the GqWo correc- 
tions to the OEPx(cLDA) starting point become k-point 
dependent, a fact already observed for GaN and II- VI 
compoundsi^ Most remarkably and in contrast to what 
we observe for GaN and InN (see Sec. IIVB 2|) the correc- 
tions are such that the band dispersion now agrees again 
with that obtained from the LDA and the LDA-f GoWo 
approach. Note also that both the band gap and 
the crystal field splitting still differ slightly between 
LDA+GoWo and OEPx(cLDA)-hGoWo for AIN [Eg-. 
LDA+GoWn : 5.95 eV, OEPx(cLDA)-HGoWo : 6.47 eV, 
Acr: LDA+GoWo : -0.252 eV, OEPx(cLDA)-hGoWo : - 
0.295 eV). Unlike for GaN and InN, experimental uncer- 
tainties do, at present, not permit a rigorous assessment. 
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which of the two GqWq calculations provides a better de- 
scription for AIN (see also Sections IIII Al and IIII Bp . 



IV. BAND DISPERSION PARAMETERS 

We will now turn our attention to band parameters 
that describe the band dispersion in the vicinity of the 
F-point: the effective masses, the Luttinger(-like) param- 
eters, and the Ep parameters. These parameters are ob- 
tained by means of the k-p -method. The k-p -method 
is a well-established approach that permits a description 
of semiconductor band structures in terms of parame- 
ters that can be accessed experimentally. Throughout 
this paper, we use four-band k-p -theory, which is typ- 
ically used to describe direct-gap materials, mostly in 
its spin-polarized form as eight-band k • p -theory. The 
k-p Hamiltonian and all relevant formulas are given in 
Appendix [X] The k - p -method is a widely accepted 
technique for, e.g., the interpretation of experimental 
data--'^-^ or modeling of semiconductor nanostructures 
and (opto-)electronic devices.— ^^i^iii^ii^iii^iii^ Its ac- 
curacy, however, depends crucially on the quality of the 
input band parameters, like effective electron masses, 
Luttinger-parameters, etc., which have to be derived ei- 
ther experimentally or from band structure calculations. 
As alluded to in the introduction, many important band 
parameters of the group-Ill nitrides GaN, InN, and AIN 
are still unknown. In particular, the band structure of 
InN is currently the subject of active research in both 
experiment and theory. 

In this paper, we use the k-p -method to derive band 
dispersion parameters from OEPx(cLDA)-|-G'oVFo band 
structures. This approach has certain advantages over 
a simple parabolic approximation around the F point. 
First, the k-p band structure is valid, not only directly 
at the F-point, but also in a certain k-range around it. 
This allows to extend the fit to larger k's and thereby 
increases the accuracy of the fitted parameters. Second, 
the k - p -method is capable of describing non-parabolic 
bands, such as the CB of InN^iiSI and can therefore also 
be applied to accurately determine values for the effective 
electron masses and Ep parameters in InN. 



A. Computational details 

For an accurate fit of the k-p parameters to the quasi- 
particle band structure a small reciprocal lattice vec- 
tor spacing is required. Since most GW implementa- 
tions evaluate the self-energy E (the perturbation op- 
erator that links the Kohn-Sham with the quasiparticle 
system) in reciprocal space, the matrix elements with re- 
spect to the Kohn-Sham wave functions (0nq|5^(erLq)l'/'nq) 
required for the quasiparticle corrections are only avail- 
able on the k-points of the underlying k-grid. A fine 
sampling of the F-point region would therefore be equiv- 
alent to using formidably large k-grids in the computa- 



TABLE V: Recommended band parameters for the 
wurtzite phases of GaN, InN, and AIN derived from the 
OEPx(cLDA)-(-G'o Wo band structures. 



par am. 


AIN 


GaN 


InN 


nic 


0.322 


0.186 


0.065 


mi 


0.329 


0.209 


0.068 


^1 


-3.991 


-5.947 


-15.803 


A2 


-0.311 


-0.528 


-0.497 


A3 


3.671 


5.414 


15.251 


^4 


-1.147 


-2.512 


-7.151 


^5 


-1.329 


-2.510 


-7.060 




-1.952 


-3.202 


-10.078 


Aj (eVA) 


0.026 


0.046 


0.175 


El (eV) 


16.972 


17.292 


8.742 


(eV) 


18.165 


16.265 


8.809 



tion. Most interpolation schemes that are frequently em- 
ployed to calculate the quasiparticle corrections for arbi- 
trary band structure points are to no avail in this case, 
because they do not add new information to the fitting 
problem at hand. Existing schemes to directly compute 
the self-energy for band structure q-points not contained 
in the k-grid (see, for instance, Ref. [Ill for an overview) 
are usually not implemented. 

In the GW space-time method^ these problems are 
easily circumvented, because the self-energy is computed 
in real-space [I]R,(r, r'; e)]. By means of Fourier interpo- 
lation, 

I]q(r,r';e) = 5]SR(r,r';e)e-*'i-^ , (1) 

R 

the self-energy operator can be calculated at arbitrary q- 
points.^'^ The matrix elements ((/'nq|5^q(e)|0nq) are then 
obtained by integration over r and r'. In this fashion 
the relevant Brillouin zone regions for the band structure 
fitting can be calculated efficiently without compromising 
accuracy. 

The k-p Hamiltonian and all parameter relations are 
given in Appendices El and |B| To determine the k - p 
Hamiltonian for a given band structure with band gap 
Eg and crystal-field splitting Aqr wc fit the parame- 
ters TOg, Ai, 7i, and Ep. This is achieved by least- 
square-root fitting of the k - p band structure to the 
OEPx(cLDA)-t-Go Wo band structure in the vicinity of F. 
For the wurtzite phases the directions S, A, T, and A 
have been included in the fit, represented by 22 equidis- 
tant k-points from F to ^ (^, ^) and 22 equidistant 
points from F to ^. For the zinc-blende phases the di- 
rections S, A, and A have been included, each with 22 
k-points from F to f (^, |). 

B. Band parameters of GaN, AIN, and InN 

The parameters obtained by fitting to the 
OEPx(cLDA)+G'oWo band structures are listed in 
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FIG. 3: Band structure wz-AlN (a), wz-GaN (b), and wz-InN (c) in the vicinity of F. The graphs show the 
OEPx(cLDA)+G'oWo band structure (black circles), the corresponding k-p band structure (black solid lines), and the k-p 
band structure using the parameters recommended by Vurgaftman and MeyepSSiiii (VM '03) (red dashed lines). 




FIG. 4: Band structure zb-AlN (a), zb-GaN (b), and zb-InN (c) in the vicinity of F. The graphs show the 
OEPx(cLDA)+GoW^o band structure (black circles), the corresponding k-p band structure (black solid lines), and the k-p 
band structure using the parameters recommended by Vurgaftman and Meyes^*^^ (VM '03) (red dashed lines). 



Tab. |V] (wz) and Tab. |Vll (zb). The resulting k-p band 
structures are plotted in Figs. [3] and |4] (black solid lines) 
together with the respective OEPx(cLDA)+G'oVKo data 
(black circles). The excellent agreement of the k-p and 
OEPx(cLDA)+G'oW^o band structures illustrates that 
the band structures of the wurtzite and zinc-blende 
phases of all three materials are accurately described by 
the k-p -method within the chosen k-ranges. Addition- 
ally, the k - p band structures based on the parameters 
recommended by Vurgaftman and Meyei^^'^ (VM 
'03) are shown (red dashed lines). As alluded to in 
the introduction, their recommendations are based on 



available experimental data and selected theoretical 
values, representing the state-of-the-art parameters 
up until the year of compilation (2003). We will also 
compare our results to more recent experimentally and 
theoretically derived parameters, (see Tab. IVIip 

In the following we will show that the parameters de- 
rived from the OEPx(cLDA)-|-GoWo calculations match 
all available experimental data to good accuracy. A com- 
parison to parameters derived by other, theoretical or 
semi-empirical, methods will be presented thereafter. 

Before we proceed, however, we would like to empha- 
size two points regarding the relation between the VB 
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TABLE VI: Recommended band parameters for the zinc 
blende phases of GaN, InN, and AIN derived from the 
OEPx(cLDA)+G'oWoband structures. The effective hole 
masses for the HH and LH band have been calculated from 
the Luttinger parameters, (see appendix [Bj 



par am. 


AIN 


GaN 


InN 


me(r) 


0.316 


0.193 


0.054 


71 






u.o 1 ^ 


72 


0.349 


0.636 


2.810 


73 


0.597 


0.977 


3.121 


'"-hh 


1.330 


0.810 


0.835 


[110] 


2.634 


1.384 


1.368 


[111] 
^hh 


3.912 


1.812 


1.738 


"hh 


0.466 


0.265 


0.080 


[110] 


0.397 


0.233 


0.078 


[111] 


0.378 


0.224 


0.077 


Ep (eV) 


23.844 


16.861 


11.373 



parameters and the effective hole masses in wurtzite 
crystals: (i) Two different sets of equations, connecting 
the effective hole masses to the Ai parameters, are used 
in the literature. Reference 1118 lists both; one is labeled 
"Near the band edge {k 0)" and the other "Far away 
from the band edge (fc is large)". The latter is widely 
used to calculate the effective hole massesJi^J^iHiii^s 
However, the experimentally relevant effective masses are 
those close to F. Thus, we use the "Near the band 
edge" equations (see Appendix |B]) throughout this work. 
Quoted values differ from the original publications in 
cases where the original work uses the "Far away from 
the band edge" equations, (ii) The Luttinger-like param- 
eters, Ai, are independent of the spin-orbit and crystal- 
field interaction parameters Aso and Acr; the effec- 
tive hole masses, however, differ for different Aso and 
Acr parameters. Only the A-band (C-band in AIN) 
hole masses can be calculated from the Luttinger-like 
parameters alone. All other hole masses depend addi- 
tionally on the choice of the spin-orbit and crystal-field 
splitting energies. ^^'^ Thus, effective B- and C-band (yl- 
and i3-band in AIN) hole masses derived from different 
sets of Luttinger-like parameters are comparable, only if 
the same Aso and Aqr values are assumed. 



1. Comparison to experimental values 

Experimentally, the band structure of a semiconduc- 
tor is accessible only indirectly, via band parameters 
like Eg, Aso, Acr, and the effective masses. Angle 
resolved direct and inverse photoemission experiments, 
which would, in principle, directly probe the quasipar- 
ticle band structure, are not accurate enough, yet, to 
determine the band structure with sufficient accuracy. 

The dispersion of the conduction band around the F 
point depends only on the effective electron masses and 
Ep parameters, which are accessible experimentally. The 



valence-band parameters, Ai, cannot be obtained directly 
experimentally, but can be related to the effective hole 
masses (see appendix [B|) . which, in turn, can be mea- 
sured. 

The available experimental values for the wurtzite 
phases are listed in Tab. lVIll For the thermodynamically 
metastable zinc-blende phases of GaN, AIN, and InN 
hardly any experimental reports on their band dispersion 
parameters are available so far. Therefore, we restrict the 
discussion to the wurtzite phases, for which experimen- 
tal data on, at least, the effective electron masses are 
available. For wz-InN also Ep has been determined, by 
fitting a simplified k-p-Hamiltonian to the experimental 
datai^2iii^ For wz-GaN, values for E p 1^5,126 g^j^^^ g^jg^ 

eral reports on the effective hole masses are availableJ^ 
Wurtzite GaN. The OEPx(cLDA)-hG'oW^o effective 

electron masses in wz-GaN {ml — 0.19 mo, = 
0.21 mo) are in very good agreement with experimen- 
tal values, which scatter around me = 0.20 too-"^ How- 
ever, our calculations predict an anisotropy of the elec- 
tron masses of about 10 %, which is larger than values 
found experimentally. (< 1% - 6%i^dS3J^) Our Ep 
values of 17.3 eV and 16.3 eV support those obtained by 
Rodina and Meyer^^s (« 18.3 eV and « 17.3 eV), rather 
than a larger value of Ep fa 19.8 eV reported recently by 
Shokhovets et alM^ 

A detailed analysis of the effective hole masses has been 
presented by Rodina et al^- Note, that only the A-band 
masses in their work have been extracted directly from 
experimental data. All other effective hole masses have 
been calculated from the A-band effective masses and 
the spin-orbit and crystal-field splitting energies within 
the quasi-cubic approximation. The effective v4-band 
masses derived in the present article {m\ = 1.88 mo and 
rrip^ = 0.33 mo) agree very well with the experimental 

values derived by Rodina et al. (mj]^ = 1.76 toq and 
= 0.35 Too). Adopting their values for the spin-orbit 
and crystal-field splitting parameters (Aso — 0.019 eV, 
Acr = 0.010 eV), we also find good agreement for the 
B- and C-band masses, (see Tab. IVII|) 

Wurtzite AIN. The available experimental data on 
the band dispersion in wz-AlN is limited to the ef- 
fective electron mass, which has been determined 
to be in the range of 0.29 to 0.45 toq.^^^ The 
OEPx(cLDA)-KGo Wo values of ml = 0.32 mo and m^ = 
0.33 mo fall within this range. 

Wurtzite InN. Experimentally derived effective elec- 
tron masses in wz-InN scatter over a wide range (see 
Tab. IVIip . The most reliable seem to be those reported 
by Wu et ai.— and Fu et al.,-^ since they explicitly ac- 
count for the high carrier concentration of their sam- 
ples and the non-parabolicity of the CB in their anal- 
ysis. Their effective electron masses of 0.05 mcr^ and 
0.07m(j21 in conjunction with values for Ep of 9.7 eV 
and 10 eV, respectively, are in good agreement with 
those derived from the OEPx(cLDA)-|-CoWo calculations 
{ml = 0.065 mo, mf = 0.068 mo and E^ = 8.7 eV, E^ = 
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TABLE VII; Band parameters of wurtzite group-Ill nitrides: Comparison to parameters from the literature. Listed are 
experimental parameters, the parameters recommended by Vurgaftman and Meyer^S, (VM '03), parameters using the LDA and 
OEPx(cLDA), parameters calculated by Carrier and Wei*'^ (CW '05) using LDA+C, and values determined using the empirical 
pseudopotential method (EPM) by several different groups. 







(recommended) 




exp. 


V iVl uo 


LDA 

(this work. 


for eomparison) 


LDA-I-C'' 


EPM 




m 


0.32 




0.29-0.45° 


0.32 


0.33 


0.38 


0.32 


0.23'', 0.24°, 0.27^ 






0.33 




0.29-0.45° 


0.30 


0.32 


0.39 


0.33 


0.24'', 0.25°, 0.18^ 


AIN 


El (oV) 


16.97 








16.41 


16.89 








Bp (oV) 


18.17 








17.45 


17.51 










3.13 






3.57 


3.46 


3.72 


3.43 


2.38'', 1.87-1. 95°, 2. 04^ 




0.69 






0.59 


0.68 


0.83 


0.68 


0.49'', 0.43-0. 48°, 0.36^ 




mil 


0.19 




0.20' 


0.20 


0.15 


0.23 


0.20 


0.14'', 0.14°, 0.16^ 




"C^ 


0.21 




0.20" 


0.20 


0.17 


0.26 


0.22 


0.15'', 0.15°, 0.12^ 




eI (eV) 


17.29 




17.8-18.7'', 19.8' 




12.93 


16.14 








-Ep (eV) 


16.27 




16.9-17.8'", 19.8' 




11.30 


14.07 






GaN 


mi 


1.88 




1.76^ 


1.89 


1.92 


2.20 


2.04 


1.89'', 2.37°", 1.45-1. 48°, 1.27-'' 




m\ 


0.33 




0.35^ 


0.26 


0.26 


0.38 


0.39 


0.26', 0.49'', 0.26-0. 27°, 0.20-'' 




m| 


0.37', 0.92" 




0.42^ 








0.85 






mi 

™c 


0.49', 0.36" 




0.51' 








0.43 






0.26', 0.19'" 




0.30' 








0.19 






0.65', 1.27" 




0.68' 








1.05 






?7ll 


0.065 


0.07 


', 0.05°, 0.04?, 0.085' 


0.07 




0.11 


0.06 


0.072", 0.10°, 0.14' 




mi 


0.068 


0.07 


', 0.05°, 0.05'', 0.085' 


0.07 




0.12 


0.07 


0.068", 0.10°,0.10' 


InN 


eI (cV) 


8.74 




10", 9.7° 






9.22 








E^ (oV) 


8.81 




10", 9.7° 






8.67 








< 


1.81 






1.56 




2.12 


2.09 


2.56-2.63^ 1.35-1.43°, 1.56^' 






0.13 






0.17 




0.21 


0.14 


0.14-0.15'", 0.18-0. 20°, 0.17^' 



"Vurgaftman and Meyer, Ref. 1261 . 
''Carrier and Wei, Ref. IstI . 
'^Silveira et aJ., Ref.ll23l. 
"^Fritsch et si., Ref . [l^. 
"^Dugdale et aJ., Ref Il24 
■fPugh et aLRef [TS^ . 
f Reference [23 and refere nces therein. 
''Rodina and Meyer, Re f.ll25l . 
'Shokhovets et ai., Ref.ll26l. 
modina et aJ., Ref [T27I. 
*Ren et aJ., Ref. [T21]. 
'Calculated using Aso = 0.019 eV 
(Ref.[l27l). 

"Calculated usineAso = 0.016 eV and Aqr 
"Wu et aJ., Ref IstI. 
°Fu et a]., Ref.lll 
^Hofmann et aJ, 



and AcR = 0.010 eV 
0.025 eV (Ref.liJ). 



'Inushima et aJ. 



Ref. 
Ref. 



12S 



"Fritsch et aJ., Ref.ll30l. 



12g 



8.8 eV). Our calculations also predict an anisotropy of the 
electron masses of about 5%. A similar anisotropy has 
been reported by Hofmann et ai.-^ (see Tab lVII[) 



2. Other •parameter sets 

Local density approximation (LDA). 
For means of comparison we have also derived band pa- 
rameters from LDA and OEPx(cLDA) calculations in the 
same way as for the OEPx(cLDA)-|-GoWo data. LDA 
band structures are frequently employed for fitting pa- 
rameter setSfiiSiii^ but we will demonstrate here, that 
the LDA is not suitable to consistently determine all pa- 
rameters for the group-III-nitrides accurately. The pa- 
rameters derived from the LDA band structures are listed 



in Tab. lVIll for the wurtzite phases of GaN and AIN. Since 
the LDA predicts InN to be metallic, no LDA band pa- 
rameters could be derived for InN. 

The effective electron masses of GaN in LDA are 
smaller than in OEPx(cLDA)+GoWo and the experi- 
ment. The effective electron masses of a given material 
are, to a first approximation, proportional to the fun- 
damental band gapj^^ Thus the underestimation of the 
effective electron masses in LDA is to some degree a nat- 
ural side effect of the underestimation of the fundamen- 
tal band gap. Additional factors (e.g. self-interaction) 
contributing to the deviation of the LDA band structure 
from the quasiparticle one were alluded to in section lll Al 

The A-band hole masses in LDA show an increased 
anisotropy; the deviation from the experimental values 
increases. 
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Despite the fact that the band gap of AIN is also sig- 
nificantly smaller in LDA than in OEPx(cLDA)-|-GoWo , 
it is still large, i.e., well above 4eV. Therefore, an effect 
on the absolute values of the effective electron masses is 
not visible, but the LDA predicts an anisotropy of the 
electron masses, with the opposite sign compared to the 
OEPx(cLDA)+GoWo calculations. 

OEPx(cLDA). 

As alluded to in section III Al band gaps in the 
OEPx(cLDA) approach open compared to LDA (cf 
Tab, imp . Following the proportionallity relationship be- 
tween the direct band gap and the conduction band ef- 
fective mass, the latter should increase in OEPx(cLDA). 
This is indeed the case, as Tab. IVIll demonstrates. They 
are, however, also larger than the conduction-band ef- 
fective masses in the OEPx(cLDA)-|-G'oWo approach, de- 
spite the fact that only in InN the OEPx(cLDA) band 
gap is larger than that in OEPx(cLDA)-|-G'oVyo . We at- 
tribute this behaviour to the approximate treatment of 
correlation in the OEPx(cLDA), which aversely affects 
the band dispersion as explained in Section IIII Dl We 
thus do not recommend the use of the OEPx or the 
OEPx(cLDA) approach alone for the determination of 
band parameters. 

LDA-plus-correction (LDA + C). 
In the LDA-I-C approacbl^I delta-function potentials arc 
added at the atomic sites, which artificially push s-like 
wave functions upwards in energy. As a consequence, the 
band gaps open, due to the admixture of cation s-states 
in the conduction band. The potentials have to be fitted 
to available experimental data, such as the fundamental 
band gaps and can be applied in an all-electron^^ii^ but 
also in a pseudopotential framework.-- Carrier and Wei 
(CW '05)^^ determined the effective electron and hole 
masses of GaN, AIN, and InN, using this method. Their 
results are also given in Tab. IVIII 

Their values for the effective electron masses of 
all three materials are in good agreement with the 
OEPx(cLDA)-|-Go Wo values. The deviations arc larger 
for the effective hole masses. This is not too surprising 
because the LDA-I-C approach predominantly affects s- 
derived bands. Since the upper valence bands around F 
are mostly of nitrogen 2p character their description will 
be closer to the LDA level, whereas the conduction bands 
feel the additional corrections. 

Empirical pseudopotential method (EPM). 
A semi-empirical way, often used to calculate band 
parameters, is the empirical pseudo potential method 
(EPM)^ii^ii24ii20 In the EPM the full atomic poten- 
tials are replaced by those of pseudo atoms, whose ad- 
justable parameters are fitted to a set of input band pa- 
rameters, typically taken from experiments. The result- 
ing band structures can then be used analogously to fit 
the parameters of a k-p Hamiltonian. Since the EPM 
depends sensitively on the input parameters, appreciable 
scatter in the reported band parameters is observed, (see 
Tab. V for a selection) 



Vurgaftman and Meyer. 
For non of the group-Ill nitrides a complete set of band 
parameters has so far been derived from experimental 
values alone. Therefore, Vurgaftman and Meyer— have 
compiled parameter sets comprising experimental and 
the most reliable theoretical values in the year 2003. 

For wz-GaN, VM'03 recommend the experimental 
value of the effective electron masses of mH = = 
0.20 mo and Luttinger-like parameters derived from EPM 
calculation by Ren et alJ^, which yield effective hole 
masses in good agreement with experimental and the 
OEPx(cLDA)-hGoWodata (see Tab. EIIl). The param- 
eter set yields a band structure that agrees well with the 
OEPx(cLDA)-fGoWoband structure for the CB and the 
two top VBs (see Fig. [3]). It deviates, however, for the 
G VB (the third valence band counted from the valence 
band maximum), where the curvatures in the EPM band 
structure are too large. 

Of all the compounds and phases discussed in this 
article wz-GaN is the best characterized experimen- 
tally. The good agreement between our quasiparti- 
cle band structures and those based on the parameter 
set recommended by VM'03 proves the quality of our 
OEPx(cLDA)-|-GoWoband structures. 

For wz-AlN the effective electron masses recommended 
by VM'03 are the averages over several theoretical values; 
the recommended VB parameters are theoretical values 
by Kim et alr^ derived from LDA calculations. These 
parameters yield a band structure, which is in good over- 
all agreement with the OEPx(cLDA)-|-GoWo band struc- 
ture (see Fig. [3^). The anisotropy of the effective electron 
masses, however, has the opposite sign, similar to our 
own LDA calculations. The similarity between VM'03 
(i.e. LDA) and OEPx(cLDA)+Go Wo in the valence band 
region (after adjusting Acr) is due to the fact that in AIN 
valence bands are shifted rigidly compared to the LDA, as 
discussed in Section fill Dl In OEPx(cLDA) alone, how- 
ever, the dispersion changes noticeably (similar to what 
was observed for InN, cf. Fig.[T]) giving rise to appreciably 
different band parameters fTab. IVTI)) . 

For wz-InN, VM'03 recommend the experimental effec- 
tive electron masses by Wu et al"^^ {ml = = 0.07 mo) 
and the EPM values from Pugh et al}^^ for the VB. The 
pseudo potentials used by Pugh et al. were designed to 
reproduce their LDA calculations, which had been "scis- 
sor corrected" to the incorrect band gap of 2.0 eV. These 
parameters are therefore to no avail from today's per- 
spective. 

V. CONCLUSION 

We have derived consistent and unbiased band parame- 
ters for the wurtzite and zinc-blende phases of GaN, AIN, 
and InN from accurate OEPx(cLDA)-|-GoWo band struc- 
ture calculations. The band parameters are in very good 
agreement with the available experimental data, proving 
the reliability of the method. We also provide reliable val- 
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ues for those parameters which have not been determined 
experimentally, such as, e.g., the band parameters of the 
zinc-blende phases of GaN, AIN, and InN or the Ep and 
VB parameters of wurtzite phases. These parameters are 
essential for understanding the physics of these materials. 
We have derived complete and consistent parameter sets 
for the description of the band structures of the group-Ill 
nitrides within k • p -theory. The k • p -method is widely 
used for modeling and simulating (opto-)electronic de- 
vices. The parameters presented in this work overcome 
the apparent lack of consistent band parameter sets for 
such simulations. 

Finally we remark that the combination of the k • p - 
with the GqWo method is not restricted to the 4x4 (8x8) 
k • p Hamiltonians discussed in this work. Since we ex- 
pect Go Wo to provide the same accuracy for the whole 



Brillouin zone, the parameters for more complex Hamil- 
tonians can be fitted in the same way. 
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APPENDIX A: K p-HAMILTONIAN 

The k-p -Hamiltonian used in the present work is based 
on the one introduced in Ref. |Tl| for zinc-blende crys- 
tals and its extens ion t o w urtzite crystals structures in 
Refs. ^ [nl [m and [Tii. It will be described in the 
following. 

Neglecting spin-orbit interaction the 8x8-Hamilton 
matrix reduces to 4x4 and can be decomposed into two 
separate matrices: 



H = Hi+H2 



(Al) 



The matrix Hi represents the pure 4x4-k-p description 
of the conduction and valence band neglecting all remote 
band contributions. For wurtzite crystals it is given by 



—iP2ky 



iP2h 

AcR + 





2mo 







A, 



CR 



hfkf_ 

2mn 






2mo ' 



(A2) 



Here, toq is the free electron mass. i?g is identical to the 
fundamental band gap for all materials with a positive 
crystal-field splitting Aqr, i.e. GaN and InN, and -t- 
|Acr| for materials with negative Aqr, i.e. AIN. The 
parameters P1/2 are proportional to the absolute value 
of the CB/VB dipole matrix elements at F. They are 
customarily expressed in terms of the Kane parameters 



E, 



Pl/2- 



P1I2 = 



2mr 



-E, 



Pl/2 



(A3) 



For zinc-blende crystals Hi simplifies through Pi — P2 
{Epi ^ Ep2) and Acr = 0. 



The matrix H2 describes the infiuences of all bands not 
considered explicitly by the 4x4- method. For wurtzite 
crystals it is defined by 



/ A'2 [kl 



A'lkl 



Ho 



Biky^kz 



L[k^ 



B2kykz 
Mik^ + M 

N!2k^k^ + N^k^ 



B2ky,K 

kyi ky 

Mikl + L'lk^ 



Bi kx kz 



ky kz 



N^k,k^ + N;,ky Msikl 



N!,k^ 



(A4) 



L'2kl j 



The parameters in H2 are defined in Ref. 11 in terms of 
optical matrix elements for zinc-blende crystals. The cor- 
responding de finit ions for wurtzite crystals can be found 
in, e.g., Ref. IllSi The parameters are related to the 



more commonly used effective electron masses, toU and 
TO;f , and Luttinger-like parameters, A\, by 

2 l^jl vaj E^ ' 
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A', = 
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mo 
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2mo 
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(A2 + A4 - A - 1) 
(Ai + A3 - 1) , 
(A2-I) , 
2A5 



2mo 

2mo 

_^ 

2mo 

_^ 

2mo 

_^ 

2mo 

i\/2A7 



\/2A6 



E, 



(A5) 



The corresponding relations for zinc-blende crystals are 
{A[ = A', =)A' 



(l; = 4 =)l' 

(Afi = M2 = A'h =) M 
{N[=N'^=)N' 



2 \m. 
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re 



(71 + 472) ■ 
(71 - 272) 

2mo 

-673 + TT 



g 
p2 



2mo 







(A6) 



Here, mo denotes the electron effective mass and 71 the 
Luttinger parameters. 

The parameters i?i/2 occur due to the lack of inversion 
symmetry in zinc-blende and wurtzite crystals. Their in- 
clusion in the k-p -Hamiltonian does not yield a noticeable 
improvement of the fit results. Therefore, they have been 
omitted throughout this work. 



APPENDIX B: EFFECTIVE HOLE MASSES 

The equations, connecting the effective hole masses to 
the Luttinger(-like) parameters are in detail 



-for wurtzite crystals — fii^ 
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For AIN the indices A, B, and C have to be interchanged: 

A^B, B^C, A. 

For zinc blende crystals it follows)^ 
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